Back

Cell Reports Medicine

Elsevier BV

Preprints posted in the last 7 days, ranked by how well they match Cell Reports Medicine's content profile, based on 153 papers previously published here. The average preprint has a 0.16% match score for this journal, so anything above that is already an above-average fit.

1
Decoding Humoral Immunity During Acute MPXV Infection via Comprehensive Serological Analysis and Antigen-agnostic Monoclonal Antibody Profiling

Zhang, Y.; Fan, J.; Wang, J.; Jiang, N.; Wan, Y.; Meng, L.; Qi, W.; Cheng, X.; Luo, K.; Zhang, T.; Li, R.; Chen, H.; Zhao, R.; Ren, Y.; Zhang, W.; Zhu, Z.

2026-08-31 public and global health 10.64898/2026.08.21.26360138 medRxiv
Top 0.5%
4.7%
Show abstract

Dissecting the complexity of antibody responses in orthopoxvirus (OPXV) infected individuals is essential for elucidating protective mechanisms and identifying candidate protective immunogens. Here, we profiled the acute humoral response in 51 mpox cases, showing distinct IgG trajectories among multiple antigens alongside the rise of plasma neutralizing activities to plateau within 6 weeks after symptom onset. Utilizing a single-cell transcriptomic and BCR sequencing based antigen-agnostic mAb isolation workflow, we further generated monoclonal antibodies (mAbs) from 254 expanded peripheral B cell clones of 3 patients. We discerned 97 specific mAbs recognizing at least 12 different OPXV proteins via integrated screening approaches, which comprised neutralizing antibodies binding unconventional viral targets and antibodies exhibiting extraordinary in vitro and in vivo anti-OPXV effects. The number of OPXV-specific mAbs recovered per donor reflected the percentage of expanded clones among circulating B cells. More interestingly, we demonstrated that the inferred unmutated common ancestors (UCAs) of neutralizing antibody clones did not necessarily react with OPXV, implying that OPXV neutralizing antibodies might frequently originate from B cells previously activated by unknown antigens. Our work establishes an efficient workflow for antigen-agnostic isolation of pathogen specific mAbs and reveals previously unclarified features of antibody responses induced by acute MPXV infection.

2
Acute Protein Responses Control SARS-CoV-2-specific Neurocognitive and General Post-Viral Sequelae

Liou, T. G.; Andrews, R. J.; Bass, B. L.; Battey, H.; Buonfiglio, L. G. V.; Cahill, B. C.; Cox, J. E.; Gibson, S.; Hartsell, S. C.; Hatton, N.; Hazel, M.; Helms, M. N.; Jensen, J. L.; Kartsonaki, C.; Kupfer, J.; Li, Y.; Lopes, F. B. T. P.; Manuel, A.; Marchetti, M.; Marvin, J. E.; Middleton, E. A.; Mimche, P.; Packer, K. A.; Paine, R.; Szczesniak, R. D.; Sturrock, A. B.; Tandar, A.; Tarbet, B.; Ulrich, A.; Warner, D.; Warren, K.; Weis, A. M.; Zimmerman, E.; Yoon, S.; Ownbey, M.; Youngquist, S. T.; Adler, F. R.

2026-08-31 infectious diseases 10.64898/2026.08.27.26361488 medRxiv
Top 0.7%
4.0%
Show abstract

Post-acute infection syndromes (PAIS) follow viral syndromes including post-acute sequelae of COVID19 (PASC) which complicates 10-25% of SARS-CoV-2 infections. These syndromes lack precise explanatory mechanisms. We studied 173 human saliva proteomes during respiratory viral syndromes, seeking associations between 44 clinically-relevant protein expression patterns and subsequent sequelae counts. Exploratory models adjusted by extensive clinical annotations found interactions between 23 acutely-responsive proteins and SARS-CoV-2 infection that inversely predicted subsequent neurocognitive sequelae. An overlapping 19 acutely-responsive proteins during any acute respiratory viral syndrome inversely predicted general fatigue-related sequelae. Altogether, 29 proteins, derived from interferon stimulated genes (ISG), were uniformly beneficial, including 13 predictive of both neurocognitive and general sequelae. The proteins suggested both shared early pathobiology and virus-specific protective responses that shaped resolution of acute disease and different PAIS. Acutely elevated protective ISG proteins associated with reduced post-viral symptoms identify investigational starting points for novel mechanisms, diagnostics and therapeutics for PASC and PAIS.

3
Multiplexed FRET-FLIM Profiling of Immune Checkpoint Interactions Predicts Response to Atezolizumab in Urothelial Carcinoma

Camacho, L.; Cacho-Navas, C.; Agüero, J.; Batmunkh, B.; Gracia, J. M.; O Sullivan, K.; Rementeria, M.; Miles, J.; Gumuzio, J.; Aguirre, F.; Martin Algarra, S.; de Andrea, C. E.; Parker, P. J.; Calleja, V.

2026-09-03 oncology 10.64898/2026.09.01.26361904 medRxiv
Top 1%
3.1%
Show abstract

Immune checkpoint inhibitors targeting the PD-1/PD-L1 axis have shown great promise in treating bladder cancer and are now part of the standard treatment for advanced disease. However, many patients still fail to respond to treatment and at present many biomarkers are assessed but have yet shown only limited results. Therefore, with the advent of combination treatments and the increase of immune related adverse event, the search for reliable predictive biomarkers is paramount. Using a multiplexed enhanced FRET-FLIM based technique (QF-Pro) we quantified the interaction of PD-1/PD-L1, CTLA-4/CD80 and TIGIT/CD155 immune checkpoints in a pre-treatment TMA of 46 patients treated with atezolizumab. The association between higher PD-1/PD-L1 ICP interaction state and treatment efficacy was demonstrated in the male sample cohort, where it identified patients with better PFS. Conversely, patients exhibiting higher CTLA-4/CD80 engagement had a worse response to atezolizumab. Remarkably, the dual assessment of patients with high PD-1/PD-L1 and low CTLA-4/CD80 allowed to identify the best responders. These results indicate that the monitoring of patients immune profile in urothelial carcinoma might be critical in identifying patients who may benefit from combination therapy.

4
Surprisal-based large language models reveal immunologic insights in lobular breast cancer

Majumder, B. P.; Linak, J. A.; Adamson, R.; Aguilera, R. L.; Agarwal, D.; Reitz, Z.; Loiselle, S.; Devarakonda, S.; Clark, P.; Paulson, K. G.; Stanton, S.

2026-08-31 oncology 10.64898/2026.08.25.26361365 medRxiv
Top 2%
2.3%
Show abstract

In large data sets discovery is often limited to pre-conceived hypotheses and data fishing. Here we tested whether systematic exploration of AI generated hypotheses could uncover clinically meaningful signals in extensively studied data. We deployed AutoDiscovery, a newly launched large language model (LLM) framework designed to search for hypotheses based on surprisal and systematically interrogate complex datasets, on The Cancer Genome Atlas breast cancer cohort. The system did not identify clinically meaningful novel findings without human input. However, a seeded warm-start run with minimal text input from an oncologist revealed multiple interesting and surprising hypotheses. Among these was that a robust immune signature was present across all subtypes of invasive lobular carcinoma (ILC) that exceeded invasive ductal carcinoma (IDC). This observation was independently validated in independent cohorts and confirmed by high-sensitivity multi-immunofluorescence tumor tissue analyses. These results suggest immunotherapy approaches should be tested in ILC including early stage ER+HER2- ILC; these patients are currently excluded from large neoadjuvant immunotherapy trials. They further demonstrate that surprisal-based hypothesis generation frameworks can extract previously unappreciated patterns from deeply interrogated cancer datasets and imply that disease domain experts working with LLMs can derive more meaningful insights from complex data than either could achieve alone.

5
A Scalable Biological Clock for Metabolic Disease Prediction from the Phenome India Cohort

Tiwari, P.; Garg, M.; Pattanayak, S.; Sarkar, I.; Roy, R.; Bhatraju, N.; Verma, A.; K, S. R.; Prakash, S.; Kumar, V. S.; Uddin, M. A.; Rawat, N.; Sahu, A.; Kumar, Y.; Leuva, P. H.; Mridha, A.; Yenamandra, V.; Singh, A. P.; Mishra, A.; Raychaudhuri, S.; Tallapaka, K. B.; Chandak, G. R.; Kulkarni, M. J.; Dharne, M.; Wahengbam, R.; Kalita, J.; Manna, P.; Subudhi, U.; Majumder, S.; Chakraborty, P.; Chaudhary, K.; Sengupta, S.; Phenome India Consortium, ; Sardana, V.; Chatterjee, S.; Ganguly, D.

2026-09-03 endocrinology 10.64898/2026.08.29.26361656 medRxiv
Top 2%
1.7%
Show abstract

Background: India has a rising incidence of chronic non-communicable diseases, making it a major healthcare burden today. Growing evidence suggests that chronic low-grade inflammation links ageing with cardiometabolic disorders, captured by the emerging concept of inflammaging. However, most evidence on biological ageing comes from Western populations, with no similar models developed for the Indian population. Given the country's distinctive genetic makeup, unique exposome, and heterogeneous NCD presentation, Western models may not capture inflammaging and its effects in the Indian population. Methods: We analysed baseline data from 4,240 adults in the Phenome India CSIR Health Cohort Knowledgebase (PI CheCK), a nationwide multi-centre cohort. Participants were stratified into eight cardiometabolic phenotype groups by BMI (Asian cut off), blood pressure and HbA1c status. We trained a Super Learner ensemble to predict chronological age in the lean normotensive-normoglycaemic reference group (n=615) using 44 plasma cytokines, sex, haemoglobin, and bioimpedance-derived visceral fat area, per cent body fat, and total body water. Performance was assessed by repeated five-fold cross-validation and in a held-out healthy test set. Calibrated biological age acceleration was then estimated in the remaining 3,625 participants. Results: Median age was 51.0 years (IQR 41.0 to 62.0) and 49.4% were female. The Super Learner outperformed elastic net and XGBoost comparators. Permutation importance identified visceral fat area, per cent body fat, CTACK, SDF1a, haemoglobin and sex as leading contributors, with body composition measures accounting for the largest share, indicating an immune-metabolic rather than cytokine-only signal. Biological age acceleration was concentrated in overweight/obese phenotypes. Lean phenotypes showed acceleration close to the reference (0.32 0.50 years). Conclusions: Cytokine and body composition measures capture a quantifiable immunometabolic ageing signal in a South Asian cohort, with acceleration driven predominantly by adiposity. External validation and longitudinal follow up are required.

6
Myelonets define spatiotemporal immunosuppressive programs in ovarian cancer

Niemiec, I.; Shabanova, A.; Ruuska, E.; Tissarinen, M.; Liang, Z.; Anandagoda, G.; Shah, S.; Kang, Z.; Junquera, A.; Salko, M.; Haltia, U.-M.; Virtanen, A.; Farkkila, A.

2026-08-31 oncology 10.64898/2026.08.26.26361128 medRxiv
Top 2%
1.7%
Show abstract

High-grade serous ovarian carcinoma (HGSC) responds poorly to immune checkpoint blockade, partly due to a macrophage-dominated immunosuppressive microenvironment. We integrated single-cell spatial proteomics and spatial transcriptomics across 50 HGSC tumors and applied SPACEstat to resolve higher-order immune communities and their transcriptional programs. We identified six immune community types, with macrophage-dominated Myelonets representing the predominant spatial pattern of immune organisation. In chemotherapy-exposed tumors, Myelonets showed coordinated lipid metabolism-immunosuppression and inflammation-MHC-II macrophage transcriptional programs, with SPP1, C1Q, VEGF, MMPs, and CCL18 linked to immunosuppressive states and fibroblasts emerging as key mediators of macrophage communication. Chemotherapy contracted large Myelonets while increasing CD8+ T-cell organization into Lymphonets. Persistent macrophage dominance within Myelonets was associated with adverse outcomes among patients who achieved a complete response to treatment. Together, we identify Myelonets as clinically relevant, multicellular immunoregulatory niches sustained by spatiotemporally coordinated macrophage programs and stromal crosstalk.

7
PCGS: biomarker and risk group identification for Pediatric Cancers via explainable Graph neural networks with Shapley values

Shi, Z.; Budhkar, A.; Amin, W.; Pollok, K. E.; Su, J.; Huang, K.

2026-09-01 health informatics 10.64898/2026.08.27.26361540 medRxiv
Top 2%
1.7%
Show abstract

Improvements in data availability, sharing, and integration, together with the development of explainable artificial intelligence (XAI) techniques, are advancing precision medicine for pediatric cancer by facilitating diagnosis, biomarker discovery, and drug development. Data sharing commons and initiatives like the Childhood Cancer Data Initiative (CCDI) provide access to pediatric-specific genomic and clinical data cohorts and improve data availability for pediatric cancer research. Based on CCDI, a scalable AI platform, Graph Artificial Intelligence for Pediatric Oncology (GAIPO), integrates various data modalities from bulk and single-cell omics data to clinical information. Such multi-modal data facilitates the training and development of advanced XAI models for pediatric cancers. We then developed an end-to-end multi-modality framework, PCGS, for pediatric cancer by incorporating omics-specific representation learning via GNN models with cross-attention fusion and multi-objective learning for downstream tasks such as classification, clustering, and survival analysis. This framework outperforms previous supervised multi-omics integration baseline approaches based on glioma and Wilms tumor cohorts and enables GNN model explainability via Shapley value-based feature attribution approaches to explain the contributions of gene-level features across various biomedical tasks, including classification and survival. Given specific background samples (e.g., age groups, sex, grades) as baselines, this explainable GNN model estimates and ranks the importance scores for input features from each omics modality. It identifies background-specific key features for biomarker discovery, risk group identification, and survival analysis in glioma and Wilms tumor, with potential applicability to other pediatric cancers.

8
Lymphodepletion mitigates anti-CAR immunity in pediatric and young adult patients with recurrent or refractory brain tumors: clinical trial results

Wang, L. D.; Oill, A. M. T.; Lindner, S. E.; Stiller, T.; Egelston, C.; Blanchard, M. S.; Mudunuri, R.; Hibbard, J. C.; Wu, M.; Sepulveda, S. M.; Peter, L.; Kilpatrick, J. L.; Stratman, J.; Mee, E. D.; Chen, D. G.; Oliveira, G.; Munoz, M.; Burmayan, A.; Wagner, J.; Dolatabadi, A. M.; Nisis, M.; Shepphird, J. K.; Sanchez, G.; Natri, H. M.; Oliver-Cervantes, C.; Feldman, L.; Aftabizadeh, M.; Arvanitis, L.; Campbell, K. M.; Cotter, J. A.; Read, J. A.; Read, J. A.; Shahani, S.; Forman, S. J.; Adam, T.; de la Nava Martin, D.; Richman, S. A.; Paul, J.; Wadden, J.; Badie, B.; Tamrazi, B.; Koschmann,

2026-09-01 oncology 10.64898/2026.08.27.26361261 medRxiv
Top 3%
1.5%
Show abstract

Outcomes for high-grade pediatric brain tumor patients remain poor, but there is optimism that chimeric antigen receptor (CAR) T cell therapy can improve prognosis. We present the results from a phase I clinical trial of IL13BBz-CAR T cells infused weekly into the cerebral ventricles in pediatric and young adult patients with recurrent or refractory brain tumors. The trial met its primary objectives of feasibility, safety, and tolerability, with one dose-limiting toxicity. 8 of 16 patients evaluable for response experienced radiographic size decreases consistent with biologic activity and with an anti-tumor response. Two patients met protocol criteria for response. Median survival for patients receiving lymphodepletion was 20.5 months from diagnosis and 6.9 months from treatment for patients with midline glioma, and 187 months from diagnosis and 7.5 months from treatment for patients with ependymoma. Importantly, patients who did not receive lymphodepletion developed anti-CAR humoral and cellular immune responses detectable in the CSF and peripheral blood, whereas patients receiving lymphodepletion had no evidence of CSF anti-CAR immunity. Taken together, these findings demonstrate the safety, tolerability, and biological activity of locoregionally-delivered IL13BBz-CAR T cells for children and young adults with CNS tumors. Moreover, we show that anti-CAR immune responses arise in patients not receiving lymphodepletion, but not in the CSF of patients receiving systemic lymphodepletion. Further investigation of adoptive cellular therapies combined with immunosuppression is warranted in this patient population. ClinicalTrials.gov registration: NCT04510051.

9
Deep phenotyping and multi-omics analyses reveal systems-wide metabolic dysregulation in a refined trisomy mouse model of Down syndrome

Saqib, M.; Chen, F.; Mistri, D. K.; Tan, L.; Wright, N.; Sarver, D. C.; Anders, R.; Aja, S.; Wong, G. W.

2026-08-29 physiology 10.64898/2026.08.26.747201 medRxiv
Top 3%
1.5%
Show abstract

Trisomy 21 or Down syndrome (DS) affects multi-organ systems across the lifespan. The presence of an extra chromosome, along with genome dosage imbalance due to triplicated genes, contributes to the DS phenotypes. Of the DS mouse models, few are aneuploid with a freely segregating extra chromosome. We previously showed that the aneuploid Ts65Dn mice exhibit metabolic deficits consistent with the metabolic profile of DS. However, the genotype-phenotype relationships in Ts65Dn mice are complicated by the presence of triplicated genes unrelated to human chromosome 21 (Hsa21). To address this issue, we leveraged a refined model, Ts66Yah, where the extra triplicated genes in Ts65Dn have been removed. Deep phenotyping and multi-omics analyses showed that Ts66Yah mice develop pronounced and widespread metabolic disturbances. Despite sexual dimorphism in weight gain, body temperature, lipid and lipoprotein profiles, hepatic injury and adipose fibrosis, both male and female Ts66Yah mice share a common phenotype of pronounced glucose intolerance and insulin resistance, reduced mitochondrial respiratory capacity in visceral fat, altered serum inflammatory cytokine profile, and dysregulated serum and liver metabolomes. Pan-tissue transcriptomes also reveal signatures of immune activation, disrupted metabolic processes and cellular respiration, altered cytokine signaling, enhanced oxidative stress, and extracellular matrix remodeling. These combined changes across tissues disrupt metabolic homeostasis more severely in Ts66Yah than in Ts65Dn mice. Several phenotypes, including glucose intolerance, insulin resistance, tissue fibrosis, and oxidative stress were further exacerbated by an obesogenic diet. This foundational data establishes Ts66Yah as a valuable reference model for the mechanistic and comparative study of metabolic dysfunction in DS.

10
Tumor-adjacent B cell infiltration stratifies recurrence risk in localized prostate cancer

Wang, B.; Mukherjee, S.; Baj, A.; Trostel, S. Y.; Lis, R. T.; Whitlock, N. C.; Ku, A. T.; Heyward, K. E.; Kartal, S.; Wang, K.; Voznesensky, O. S.; Calagua, C.; Siddiqui, J.; Martin, R. S.; Kollath, L. A.; Custer, J.; Michael, P. D.; Kunju, L. P.; Lake, R.; Harris, C. C.; Aldape, K. D.; True, L. D.; Tatsuoka, C.; Fertig, E. J.; Chinnaiyan, A.; Gurram, S.; Pinto, P. A.; Weiner, A. B.; Morrissey, C.; Salami, S. S.; Einstein, D. J.; Balk, S. P.; Sowalsky, A. G.; Ruppin, E.

2026-08-31 oncology 10.64898/2026.08.29.26361718 medRxiv
Top 3%
1.5%
Show abstract

Background: Biochemical recurrence (BCR) occurs in 20-40% of men after radical prostatectomy. Existing postoperative recurrence risk tools based on PSA and pathology are clinically useful but show only moderate and variable discrimination, highlighting the need for biomarkers that improve risk stratification and consequent treatment decisions. We hypothesized that the prostate microenvironment, including both the tumor and non-cancerous adjacent tissue, may contain prognostic features associated with adverse postoperative PSA outcomes. Methods: We assembled a cohort of matched tumor-adjacent benign and tumor prostate tissue from 243 men across three institutions to establish a discovery cohort (n=123; 43 postoperative PSA events, 35%) and validation cohort (n=120; 46 events, 38%). For primary binary analyses, a postoperative PSA event included BCR, defined as two consecutive postoperative PSA values >=0.2 ng/mL, or PSA persistence. We performed RNA sequencing of matched tumor-adjacent benign and tumor tissues, quantified immune signatures, and developed an integrated model combining the adjacent-tissue B-cell signature, preoperative PSA, and radical prostatectomy Gleason score (BRIGADE). CAPRA-S-adjusted Cox analyses excluding recurrence-time-0 cases evaluated time to BCR, and CD19 multiplex immunofluorescence provided tissue-level confirmation (n=10). Results: In prostatectomy specimens, tumors from patients without a postoperative PSA event were enriched for B-cell transcriptional programs, whereas tumors from event-positive patients showed elevated proliferation signatures. B-cell-related transcriptional programs were correlated between tumor and adjacent tissue. Tumor-adjacent benign B-cell scores were higher in no-event cases and discriminated postoperative PSA-event status in PCBN discovery (AUC 0.63) and BM validation (AUC 0.81) cohorts, outperforming numerous other immune-related signatures. In CAPRA-S-adjusted Cox sensitivity analyses excluding recurrence-time-0 cases, higher adjacent-tissue B-cell activity was associated with reduced recurrence risk in PCBN (HR 0.42, 95% CI 0.19-0.94; BH-adjusted p=0.035) and BM (HR 0.54, 95% CI 0.30-0.95; BH-adjusted p=0.034). Tissue-based validation showed that CD19+ B-cell density in adjacent benign tissue was higher in no-event than event-positive patients (median 0.1145 vs 0.0471; p=0.008). BRIGADE achieved an AUC of 0.68 in cross-validation and 0.83 in independent validation, compared to AUCs of 0.54-0.63 and 0.44-0.78 for the tested clinical predictors, respectively. At the fixed classification threshold, the validation-cohort odds ratio for BRIGADE was 2.75. The adjacent B-cell score remained associated with lower odds of a postoperative PSA event after adjustment for PSA and Gleason score. Conclusions: B-cell infiltration in tumor-adjacent benign prostate tissue may complement existing clinicopathologic models for stratifying adverse postoperative PSA outcomes and subsequent BCR after radical prostatectomy. The transcriptomic signal was recapitulated by CD19-based tissue staining, supporting further development of a pathology-based assay.

11
The Role of Distress-related Metabolic Dysfunction in Ovarian Cancer Development: a pooled case-control study

Lin, N.; Balasubramanian, R.; Menichetti, G.; Eliassen, H.; Trabert, B.; Avila-Pacheco, J.; Townsend, M. K.; Terry, K. L.; Clish, C. B.; Tworoger, S. S.; Zeleznik, O. A.

2026-08-31 epidemiology 10.64898/2026.08.27.26361473 medRxiv
Top 3%
1.2%
Show abstract

Background: Evidence suggests chronic distress influences ovarian cancer (OC) etiology and metabolomic profiles. Here, we evaluated the association of a metabolite-based distress score (MDS) and OC risk. Methods: We included two matched case-control studies nested within the Nurses' Health Studies (N=584) and the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial (N=348). Metabolites were measured 3-27 years before diagnosis using liquid-chromatography tandem mass spectrometry. We examined the association of quintiles of MDS and 19 constituent metabolites with OC risk using unconditional logistic regression and stratified by tumor histotype, menopausal status, and age at diagnosis. Results: We observed women in the highest versus lowest quintile of MDS had an increased OC risk (OR=1.62,95%CI=1.03-2.54,ptrend=0.07), and type 2 tumors (OR=1.71,95%CI=1.03-2.83,ptrend=0.11). Associations were suggestively stronger for premenopausal and <69-year-old women, and driven by pseudouridine, and N2,N2-dimethylguanosine. Conclusion: Our findings suggest chronic distress-associated metabolic dysregulation may represent a novel OC risk factor, especially among younger women.

12
Identification of genetic variants in Pfs25 and functional evaluation in mosquito infection

Orfano, A.; Cisse, A.; Guo, Y.; Han, L.; Fikadu, N.; Thiam, L. G.; Ba, A.; Li, R.; Pouye, M. N.; Mangou, K.; Moore, A. J.; Sene, S. D.; Diallo, F.; Ngom, E. M.; Sadio, B.; Mbengue, A.; Membi, C.; Ngasala, B.; Bazie, T.; Some, F. A.; Olson, N.; Patel, S. D.; Shapiro, L.; Parikh, S.; Foy, B. D.; Cappello, M.; Vigan-Womas, I.; Premji, Z.; Dabire, R. K.; Ouedraogo, J.-B.; Sheng, Z.; Bei, A. K.

2026-08-31 infectious diseases 10.64898/2026.08.25.26361130 medRxiv
Top 4%
1.1%
Show abstract

Transmission-blocking vaccines (TBVs) are a promising strategy to reduce malaria transmission by targeting parasite stages within the mosquito. However, parasite genetic diversity may limit vaccine efficacy. We used next-generation amplicon deep sequencing to identify non-synonymous single nucleotide polymorphisms (SNPs) in Pfs25 from 184 Plasmodium falciparum isolates from Senegal, Tanzania, Ghana, and Burkina Faso. Prioritized SNPs were introduced into P. falciparum via CRISPR-Cas9. For the G116C variant, gametocyte development was evaluated by microscopy and qPCR, and mosquito infectivity was assessed by SMFAs. We identified 26 SNPs, including 24 novel variants. Functional assays showed that the Pfs25 G116C mutation did not affect gametocyte development or exflagellation. SMFA showed no significant differences in oocyst prevalence or intensity between mutant and WT parasites. These findings highlight the importance of integrating genetic surveillance with functional validation to guide the development of effective transmission blocking interventions

13
Hybrid risk scores integrating polygenic and clinical variables for endometriosis prediction

Goroshchuk, O.; Koller, D.

2026-09-03 epidemiology 10.64898/2026.08.31.26361798 medRxiv
Top 4%
1.0%
Show abstract

Background: Endometriosis affects approximately 10% of reproductive-age women and is associated with substantial diagnostic delay and heterogeneous symptom presentation. Prior machine-learning prediction models have relied on comorbidity data alone or on small candidate-variant genetic scores, with inconsistent or incompletely reported performance. No study has combined a well-powered, multi-ancestry polygenic risk score (PRS) with environmental, reproductive, and symptom data in a single hybrid model. We developed and evaluated hybrid risk-prediction models integrating a genome-wide, multi-ancestry PRS with clinical and symptom data for endometriosis in the US-based All of Us Research Program. Methods: Among 69,376 participants (15,382 endometriosis cases, 53,994 controls) across six genetically inferred ancestry groups, we computed individual-level PRS values using PRS-CS weights derived from an independent, multi-ancestry GWAS. Five nested logistic regression, random forest, and XGBoost models progressively added age, ancestry, and within-ancestry genetic principal components (Model 1), environmental and reproductive factors (Model 2), symptom and comorbidity indicators (Model 3), all covariates combined (Model 4), and PRS x environment interactions (Model 5). Performance was assessed by AUROC in a held-out test set and 5-fold cross-validation, with class-weighted, Youden-optimized thresholds used for sensitivity, specificity, and predictive values; permutation importance identified top contributors. Pairwise AUROC differences were tested with a Holm-corrected DeLong-type test. Results: Discrimination improved from AUROC 0.63 (PRS, age, ancestry, principal components) to 0.72 for the full model, driven mainly by symptom and comorbidity data. XGBoost consistently outperformed logistic regression and random forest. The PRS ranked among the top individual predictors by permutation importance in nearly every model, alongside age, while genetic and demographic information alone gave only modest discrimination, and PRS x environment interactions did not improve on environmental factors alone. Threshold optimization yielded balanced sensitivity and specificity (~0.67/0.65) versus near-zero sensitivity at a default threshold. Conclusions: Combining the PRS with symptom and comorbidity data gave the best discrimination compared to solely a well-powered, multi-ancestry PRS as a predictor of endometriosis. This study clarifies both the promise and current limits of hybrid genetic-clinical prediction for endometriosis and points to symptom-based phenotyping, molecular subtyping, and external validation as priorities.

14
AURORA: Analysing and understanding responses to oncological regimens with artificial intelligence

Lebmeier, A.; Lindner, T.; Karl, C.; Schöler, T.; Rank, A.

2026-09-02 health informatics 10.64898/2026.08.30.26361778 medRxiv
Top 4%
1.0%
Show abstract

Background: Immunochemotherapy (ICT) is considered standard in regards to care for small-cell lung cancer (SCLC) in extensive stages, yet reliable biomarkers for treatment response remain elusive. While previous univariate analyses suggest specific peripheral lymphocyte subsets correlate with survival, the systemic immune response involves complex, multivariate interactions that require advanced analytical approaches. Methods: This paper analysed high-dimensional flow cytometry data from 32 patients with stage IV SCLC treated with carboplatin, etoposide, and atezolizumab. Peripheral blood was analysed at baseline (V0) and longitudinally during treatment. To identify potential early predictive biomarkers and mitigate sample attrition in later cycles, we focused on baseline and measurements after two cycles of ICT (V1). We employed a rigorous machine learning framework utilising nested cross-validation, bootstrapping, and permutation-based statistical testing to evaluate eleven different regression and survival models. Results: Under model-appropriate metrics, regressors did not generalise (R2 <0); conversely, censoring-aware Random Survival Forests (RSF) successfully extracted robust prognostic signatures. Baseline immune profiles (V0) achieved a concordance index (C-index) of 0.66 (p= 0.015), while dynamic changes from V0 to V1 ({triangleup}V) achieved a C-index of 0.65 (p= 0.022). Crucially, absolute values measured after two cycles of ICT (V1) yielded no significant signal (p= 0.445). Feature importance analysis confirmed the prognostic value of Th17 normalisation and identified Naive Regulatory T cells and Memory B cells as candidate components. Conclusion: Machine learning validation confirms a predictive signal in the peripheral immune profile of SCLC patients. Early dynamic shifts in the balance between regulatory and effector immune arms are associated with prognosis, contrasting with the lack of signal in absolute counts after two cycles of ICT. These findings establish a proof of concept for multivariate liquid biopsy immune profiling, warranting confirmation in larger cohorts and highlighting the necessity of integrating systemic and tumour-intrinsic data.

15
Diversification without convergence: national childhood respiratory pathogen spectra diverge as they diversify, 1990-2023

Li, D.; Feng, Q.; Zhang, Y.; Chen, H.; Wang, X.; Shen, C.

2026-09-03 pediatrics 10.64898/2026.09.01.26361890 medRxiv
Top 4%
1.0%
Show abstract

Background National childhood respiratory pathogen spectra are diversifying nearly everywhere - within-country diversity rose in 203 of 204 countries between 1990 and 2023 - yet whether countries are diversifying toward a common spectrum or along divergent paths is unknown. We quantified between-country compositional distance of national pathogen spectra over the same period. Methods We built national pathogen share vectors from Global Burden of Disease Study 2023 lower respiratory infection etiologic attributions (26 pathogens, 204 countries, ages 0-19 years) at five timepoints spanning 1990-2023. Between-country distance was measured as all pairwise Jensen-Shannon divergences (JSD; primary) and Bray-Curtis dissimilarities, with Baselga and Jaccard decompositions; robustness was assessed across metrics, pathogen panels, low-count thresholds and a balanced panel of 107 countries. Results Mean pairwise JSD rose from 0.0084 in 1990 to 0.0283 in 2023 (+238%; trend p = 0.030), peaking in 2021 (+283%) with a partial 2023 pullback. Bray-Curtis dissimilarity rose +120% and the balanced panel +423%. Divergence was entirely balanced variation (share reallocation), with spectrum richness rising from 18.5 to 21.1 of 26 pathogens. Dispersion rose fastest for influenza (coefficient of variation 0.03 to 0.55) and respiratory syncytial virus (0.08 to 0.48). Within-region distance rose in every computable GBD super-region (five of seven): divergence occurs within regions, not between blocs. Conclusions National spectra are re-sorting along country-specific axes as vaccine-preventable dominance recedes at different speeds. Diversification is universal, but convergence is absent: the transition at the etiologic-spectrum level is asynchronous and path-dependent, with implications for empirical treatment policy and pathogen surveillance.

16
Antibody profiles across H5N1 and previously circulating viruses are highly dynamic and age- and imprint- independent

Beukema, M.; Vermeulen, E.; de Vries-Idema, J.; Huckriede, A.; Joshi, M.

2026-08-31 infectious diseases 10.64898/2026.08.26.26361396 medRxiv
Top 4%
1.0%
Show abstract

The increasing incidence of H5N1 influenza virus transmission from animal species to humans has heightened concerns about an imminent H5N1 pandemic. Prior studies using recombinant hemagglutinin and neuraminidase proteins have reported age-dependent cross-reactivity to H5N1, attributed to immune imprinting from an individual's first influenza virus exposure. However, whether this pattern holds when using whole inactivated virus (WIV), capturing antibodies against diverse viral proteins, and is stable over time remains unknown. We therefore aimed to determine whether H5N1 cross-reactivity of pre-existing antibodies to whole virus follows an age-dependent or imprinting-specific pattern, and whether this pattern is stable over a five-year period. To this end, we measured serum antibody levels in adolescents, adults and seniors by ELISA using whole inactivated H5N1 virus as antigen rather than purified proteins. Detectable, albeit generally low, levels of H5N1-reactive antibodies were present in most individuals, irrespective of age. Comparison of antibody levels against H5N1 with those to five historical influenza virus strains revealed a consistent positive correlation between H5N1-reactive antibodies and responses to the H1N1pdm09 strain A/California/7/2009 (CA), across all age groups. Using unbiased clustering of antibody titers against H5N1, CA, and the H3N2 strain A/Perth/16/2009 (PE), we identified seven distinct age-transcending antibody profiles. These profiles covered individuals with varying titers to all three included viruses but also identified individuals with high anti-CA levels, yet low anti-H5N1 levels and vice versa. Moreover, despite stable antibody levels over a five-year interval in the study population, individual antibody levels and profiles fluctuated considerably over this period. Taken together, our results confirm the presence of H5N1-reactive antibodies in human sera and their association with previously circulating strains. However, they also caution against inferring antibody levels against a new strain based solely on responses to antigenically related strains and highlight the limitations of extrapolating immune status from single timepoint measurements.

17
Both ageing and frailty status impact vaccine-induced transcriptomic profiles and subsequent humoral immunity: results from the VITAL cohort

Joshi, M.; Carre, C.; Cevirgel, A.; Bijvank, E.; Chabaud-Riou, M.; Courtois, V.; Chautard, E.; Larocque, D.; Burny, W.; Beckers, L.; Buisman, A.-M.; Rots, N.; van der Heiden, M.; van Beek, J.; van Sleen, Y.; van Baarle, D.

2026-08-31 allergy and immunology 10.64898/2026.08.26.26361408 medRxiv
Top 4%
1.0%
Show abstract

Vaccine responses vary across individuals due to differences in ageing and health status. Using transcriptomic profiling, we analyzed early gene expression profiles after influenza (QIV) followed by pneumococcal (PCV13) vaccination in 148 participants spanning young, middle-aged, and older adults. The two vaccines induced distinct immune signatures: QIV elicited innate and interferon immune activation, while PCV13 triggered inflammation-based responses. Older adults showed weaker but similar transcriptomic profiles compared to young adults. Among older adults, frailty, in addition to age, was strongly associated with reduced innate responses. In addition, we identified associations between early-stage transcriptomic profiles and later-stage antibody responses for QIV; however, no such associations were observed for PCV13. Importantly, observed group differences arose not from altered immune modules but from differences in the magnitude of gene expression, paving the way for immune-boosting interventions to enhance early gene expression in at-risk populations.

18
Impaired memory B-cell formation after mRNA-based COVID-19 booster vaccination in patients with inflammatory bowel disease receiving anti-TNF treatment

Gill, P. A.; Bradbury, L. R.; Wang, A.; Hogg, J.; Demase, K.; McKenzie, J.; Fryer, H. A.; Geers, D.; Zaeck, L. M.; Boo, I.; Hogarth, M. P.; Drummer, H. E.; de Vries, R. D.; O'Hehir, R. E.; Sparrow, M. P.; van Zelm, M. C.

2026-09-02 allergy and immunology 10.64898/2026.08.28.26359302 medRxiv
Top 5%
0.8%
Show abstract

Background: Patients receiving anti-TNF treatment for chronic inflammatory disease display impaired antibody responses, but it remains unclear how immune memory formation is affected. We evaluated antibody responses and memory B cells (Bmem) after COVID-19 booster vaccination in inflammatory bowel disease (IBD) patients receiving anti-TNF treatment. Methodology: Blood was sampled at baseline, 1, and 6 months after WH1/BA.5 bivalent or XBB.1.5 monovalent vaccination from 27 IBD patients receiving intravenous anti-TNF and 44 controls. Neutralizing antibodies were measured using an infectious virus assay. SARS-CoV-2 spike receptor binding domain (RBD)-specific serum IgG was quantified by ELISA, and RBD-specific Bmem were immunophenotyped by flow cytometry using recombinant proteins from ancestral, Omicron BA.1, BA.5, XBB.1.5, and JN.1 variants. Results: Serum IgG to vaccine RBD and neutralizing antibodies in patients increased pre to 1 month post-vaccination, but were lower than controls. Ancestral-, BA.5- and XBB.1.5-specific Bmem increased after vaccination but were significantly lower in patients than controls. Within RBD-specific Bmem, frequencies of recently activated CD21lo cells were increased after vaccination, and were higher in patients than controls. Fewer antigen-specific Bmem in patients expressed IgG4, and more expressed IgG3 or IgD following vaccination. Following vaccination, more RBD-specific Bmem recognized multiple viral variants. However, patients had fewer Bmem that could bind to subvariants than controls. Conclusion: Antibody and Bmem responses to COVID-19 booster vaccination in anti-TNF-treated IBD patients displayed reduced capacity, durability and cross-reactivity, suggesting impaired immune memory for protection against breakthrough infection. This supports the recommendation for annual booster vaccination to prevent severe disease and viral spread.

19
Prospective In-silico Simulation of the VESALIUS-CV Trial Using Biomedical Knowledge Graph and Real-World Data-Driven AI Modeling

Perlman, A.; Goldstein, N.; Goldman, M.; Shapiro, M.; Barash, E.; Bar, A.; Raveh, T.; Tordjman, E.; Schussheim, H.; Dormont, F.; Matalon, O.

2026-08-31 cardiovascular medicine 10.64898/2026.08.26.26361436 medRxiv
Top 5%
0.8%
Show abstract

Background. Cardiovascular-outcomes trials are lengthy, costly, and associated with substantial uncertainty prior to readout. In-silico trial simulation using real-world data (RWD) has emerged as a potential tool to support earlier decision-making; however, evidence of prospective predictive validity, generated prior to trial result disclosure, remains limited. Methods. We applied a semi-mechanistic machine learning framework integrating real-world patient data with biologically informed drug representations to prospectively simulate the VESALIUS-CV trial evaluating evolocumab versus placebo. The simulation model was trained on a combination of patient-level real-world data and a drug-centric knowledge graph and validated for both patient-level and trial-level retrospective predictive performance. The model was then used to simulate VESALIUS-CV before public disclosure of trial results, using a locked model and prespecified eligibility criteria and primary endpoint aligned with the clinical protocol. A patient-level time-to-event model was used to generate virtual trial arms, from which cumulative incidence curves, hazard ratios, confidence intervals, and p-values for major adverse cardiovascular events (MACE) were estimated. Results. In retrospective validation, the model demonstrated strong patient-level discrimination, with time-dependent ROC-AUC values ranging from 0.80 to 0.90 across follow-up horizons. For trial-level validation, 22 randomized cardiovascular-outcomes trials were simulated, and hazard ratios for 3-point MACE across 24 between-arm comparisons showed consistent directional agreement and quantitative correlation with published results such that the model accurately predicted trial success, achieving an F1 score of 0.83, with precision of 0.79 and sensitivity of 0.89. In a fully prospective application, the simulation predicted a statistically significant reduction in 3-point MACE with evolocumab versus placebo, estimating a hazard ratio of 0.78 (95% CI, 0.70-0.87) at 54 months. These predictions were consistent with the subsequently reported VESALIUS-CV results, which demonstrated a hazard ratio of 0.75 (95% CI, 0.65-0.86) at 55 months of median follow-up. Conclusions. In a fully prospective setting, a RWD-driven, AI-based simulation accurately predicted the direction, magnitude, and temporal dynamics of treatment effects observed in the VESALIUS-CV trial. These results demonstrate that in-silico trial simulation can anticipate clinical outcomes in the prospective setting, supporting its use as a complementary tool for early decision-making, trial design optimization, and de-risking in cardiovascular drug development.

20
Genome Profiling of Actionable Cancer Targets (NYU LG-PACT) for Clinical Patient Molecular Diagnostics and Treatment

Yang, Y.; Vasudevaraja, V.; Serrano, J.; Mohamed, H.; Kelly, S.; Jour, G.; Gindin, T.; Park, K.; Jones, D.; Feng, X.; Pinnell, J.; Mclennan, S.; Tin, M. Y.; Tsirigos, A.; Snuderl, M.; Wrzeszczynski, K. O.

2026-09-01 oncology 10.64898/2026.08.27.26361341 medRxiv
Top 5%
0.6%
Show abstract

Next-generation sequencing (NGS) for the detection of somatic variants has become the method of choice in a variety of molecular oncology fields and in the clinic. Its use ranges from sequencing entire tumor genomes and transcriptomes to targeted clinical diagnostic gene panels. The NYU Langone Genome PACT (Profiling of Actionable Cancer Targets, LG-PACT) assay is a qualitative in vitro diagnostic test that uses targeted next generation sequencing (NGS) of formalin-fixed paraffin-embedded (FFPE) tumor tissue matched with normal specimens from patients to detect gene alterations in a targeted panel covering 606 genes and the TERT promoter. Indications for testing are cancer (solid tumors and hematological malignancies) where a mutational profile from multiple genes would be informative for disease stratification, prognosis, or treatment options including targeted therapies and eligibility for clinical trials. The test is intended to provide information on somatic mutations including point mutations, small insertions/deletions (indels), and copy number aberrations for diagnostic and treatment decisions. LG-PACT is a United States Food and Drug Administration (FDA) cleared diagnostic test (510K: K202304). The clinical interpretation of sequencing data of molecular tumor markers from NGS encompasses automated variant calling tools with human interpretation. This final mostly manual review of data step is intensive, involving highly trained scientists, encompassing literature review, interpretation and clinical tier classification by pathologists, who then provide a complete molecular diagnostic report to the treating oncologists. We provide analysis of 1339 clinical genomic profiles from 31 different cancers and their subtypes, comprising of central nervous system (CNS) 792 (59%) cases (incl. meningioma, glioma and glioblastoma), with 267 (20%) cases predominantly of lung, pancreatic and colorectal and 280 of others (21%). Here, we present the technical challenges of validating an NGS oncological diagnostic targeted assay for clinical grade accuracy and sensitivity for patient care. We show how copy number alterations provide a more comprehensive description of the tumors genomic profile. We then outline the utility of targeted panel sequencing based on certified pathologist selection of reportable variants for our current patient cohort. Where analysis of variant detection has led to 49.4% (661/1339) of our clinical tumor samples containing mutations in known therapy targeted genes, 35.6% (477/1339) with mutation detected in other genes, and 15% (201/1339) cases being negative.